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Abstract. The law of reflection and Snell's law are among the tenets of geometrical 
optics. Corrections to these laws in wave optics are respectively known as the angular 
Goos-Hanchen shift and Fresnel filtering. In this paper we give a positive answer to the 
question of whether the two effects are common in nature and we study both effects in 
the more general context of optical beam shifts. We find that both effects are caused 
by the same principle, but have been defined differently. We identify and discuss the 
similarities and differences that arise from the different definitions. 
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1. Introduction 

It is not uncommon in physics that the same effect is discovered independently by 
different people and, as a consequence, may be known under two or more different 
names owing to differences in language and scientific tradition. A prominent and 
relevant example for this article is Snell's law, which in French speaking countries is 
more commonly know as Descartes 's law, to remind us that it was Rene Descartes who 
first published it in modern times [1] (although Ibn Sahl, Willebrord Snellius and others 
derived the result before him [2]). 

The deviation from the predictions of Snell's law in wave optics is known as 
'Fresnel filtering' (FF) [3j and describes a correction to the angle of refraction, which 
is particularly prominent for focussed beams incident close to the critical angle of total 
internal reflection. The origin of this effect is easily understood from the following 
simple picture (see Fig [T]): assuming a symmetric angular spectrum and that the 
central wave vector of the focussed beam be incident exactly at the critical angle, all 
plane waves of the incident beam with an angle of incidence larger than the critical 
angle are fully reflected. On the other hand, all plane waves incident below critical 
incidence are transmitted. The magnitude of the effect is then the difference between 
the angle of refraction expected from geometrical or ray optics and the far field angle 
of the transmitted field. If the incident angular spectrum is large, so is the transmitted 
spectrum, which explains intuitively why the effect of Fresnel filtering depends on the 
angular width of the incident beam. 

Independently from Fresnel filtering a family of optical beam shift effects has 
been studied mainly on reflection [4]. In contrast to Fresnel filtering, which has been 
discovered in the context of two dimensional, scalar fields in optical microcavities [51 E], 
optical beam shifts have been explored in all three dimensions. Confined within the 
plane of incidence and therefore comparable to the 2D fields in microcavities, are the 
spatial and angular Goos-Hanchen shift (GH) [jj |8] and it is the latter which concerns 
us here as a related effect to FF. In its simplest form the angular GH shift pertains to 
generic incidence angles, that is not for the critical angle or the Brewster angle, although 
these special cases have been studied too [9j HOj [H] . Using the same simple picture as 
above, but for generic incidence (see Fig. [2]), the angular GH shift can be explained as 
a weighting in the Fourier spectrum of the incident beam upon reflection, which leads 
to a shift of the mean angle in the reflected beam. Unsurprisingly, the magnitude of the 
angular GH shift also depends on the angular spectrum. 

One difference between the two effects is based on the different definitions as 
peak and mean angle of the transmitted field. Another distinction between FF 
and the angular GH shift is rooted in their importance for the different scientific 
communities. For example, the presence of regular, triangular-type resonator modes 
in chaotic microcavities of spiral shape could only be explained by the FF effect at near- 
critical reflection [12] . Such quasi-scarred resonances [HI [15] had been observed in wave 
calculations and, implementing the FF effect into an extended ray dynamics, could now 
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Figure 1. (Colour online) Schematic of Fresnel filtering at critical incidence. A 
light beam is incident from the region of higher refractive. The central wave vector 
impinges on the interface precisely at the critical angle of total internal reflection. 
All plane waves of the incident angular spectrum (green gradient) with larger angles 
of incidence (crossed) are fully reflected, whereas plane waves with smaller angles of 
incidence are transmitted (blue gradient). The gradient indicates the amplitude of the 
angular spectrum. Right at the critical angle the transmittance is zero which is why 
the transmitted beam has the maximum shifted away from grazing transmission. For 
clarity the reflected beam is shown as arrows only. 

be understood in terms of ray optics. (Notice that implementing the FF effect alters 
the ray optics from Hamiltonian to become non-Hamiltonian [T21 [13], an interesting 
feature that is, however, beyond the scope of the present paper). This example of FF in 
quasi-scarred modes illustrates nicely two properties of the effect: (i) The FF effect is 
noticeable especially around the critical angle, which, given a refractive index of n = 2 
inside the cavity, happened to be the inner angle of the triangular-type quasi-scarred 
modes, (ii) The FF effect depends on the resonance wave number, implying a slight 
rotation between adjacent resonances in order to adjust the resonance geometry to the 
slight change in the FF correction. 

In beam optics, on the other hand, the angular shift had been found to occur in 
the same setting as the spatial GH shift, but for incidence below the critical angle 
fT6l IT7j . As the angular GH shift is a deflection of the beam's propagation direction 




Figure 2. (Colour online) Schematic of the angular Goos-Hanchen (GH) shift. 
A lightbeam is reflected from a region of higher refractive index. Upon reflection 
each plane wave acquires a real weighting by virtue of the angle-dependent reflection 
coefficient. This changes the angular spectrum and as a consequence the mean angle 
of the reflected beam is different from the angle of incidence. 



the displacement of the beam caused by this effect scales with the distance from the 
interface. This is why in situations where both spatial and angular GH shifts occur, for 
example in metallic reflection [T8] . the contribution of the angular shift to the overall 
displacement will be become dominant when looking sufficiently far down the beam 
[T9] . Recently, it has been found that for higher order Laguerre- Gaussian beams, or 
more generally any vortex beam, the total spatial shift can be explained by a mixing of 
the spatial and angular shifts of the fundamental beam [20l ETJ [22] . 

FF and the angular GH shift are both counterparts to spatial GH shift. Either as 
two independent directions in phase phase [23], or as real and imaginary part of the 
same complex quantity namely the logarithmic derivative of the reflection coefficients 
f24l |9j [8j [25] . As such the formulas for the angular and spatial GH shift on reflection 
are 

1 



A a GH = a 2 Re 



and 



-Im — 

k r 



(i) 
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where r is the reflection coefficient [271 [28] and the prime denotes the derivative with 
respect to the incident angle at the central angle 9 . These shifts have their natural 
units, given by the angular width or variance a 2 for the angular shift and the inverse of 
the wavenumber k for the spatial shift. 

In this paper we focus on the differences and similarities between the angular Goos- 
Hanchen (GH) shift and Fresnel filtering. The arrangement of the paper follows our 
logical outline of the argument. In the first part of the next section we derive an explicit 
expression for the angular GH shift on transmission for generic incidence which we 
compare to the known formula for the corresponding shift on reflection. In the second 
part we find an explicit expression for critical incidence and we compare it with the 
known results from FF at critical incidence. Section [3] contains derivations of analytical 
formulas for FF at a generic angle of incidence, which we use for a concluding comparison 
in a discussion in section [U 

2. Angular GH on transmission 

Our derivation for the expression of the angular GH shift in transmission relies on the 
concept of the virtual beam as introduced in [25J. However, whereas for the case of 
reflection the virtual beam is obtained by a specular reflection of every plane wave, 
that is the reflection coefficient is r = ±1 throughout, for transmission the virtual 
refracted beam is obtained by applying Snell's law for every plane wave according to 
the local angle of incidence. The spatial and angular shifts on transmission are then 
the differences to the virtual reflected beam due to the inclusion of the appropriate 
transmission coefficient. 

In 2D we can treat the polarisation orthogonal to plane of incidence (s or TM) 
and parallel to it (p or TE) separatelyjf] and the incident and virtual transmitted beam 
are the same for both polarisations [27] ; the only place where the polarisation enters 
is the choice of the appropriate transmission coefficient. We adopt the notation of a 
reduced refractive index n = n^jn^ where n\ is the index of the medium in which the 
incident beam and reflected beam propagate. Using x for the spatial coordinate along 
the interface and z for the coordinate normal to it, we parametrize the incident beam 
by means of the angle of incidence 9 for each plane within the spectrum a (9) 



The virtual transmitted beam is constructed by changing the direction of every plane 
wave in accordance with Snell's law sin 9 = n sin r while maintaining the parametrization 



X As there is always a confusion between the different labels for polarisation, we state explicitly that p 
(TE) is the polarisation for which there exists a Brewster angle in partial reflection. Our identification 
of s being TM and p being TE follows common usage in the microcavity community where the reference 
plane is the 2D resonator and not the plane of incidence. 




(2) 



of the beam by 9 (see Fig. [3]): 




(3) 
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If the incident beam is in the optically thicker medium, the virtual beam includes 
evanescent waves. For the calculation of the angular GH shift, however, we only include 
propagating plane waves, which can be enforced by setting the upper bound of the 
integration to the critical angle 9 C . The real transmitted beam is obtained from the 
virtual beam by inserting the appropriate transmission coefficients for s (TM) and p 
(TE) polarisation [28] 

t »= 2 "»» _ (TM), (4) 

cos 9 + v n 2 — sin 9 

2ncos# /rTmX , v 

t P = ^= (TE), (5) 

n 2 cos 9 + v n 2 — sin 2 9 



^ = / d9a(9)t a (9)e ik ( xsiQd - zVn2 - sin2d ) , (6) 



which yields 



where a = s(TM),p(TE). 

Following the treatment in [25] we calculate the mean angle of the transmitted beam 
ifj t as centroid of the transmitted angular spectrum. Unlike for reflection, however, we 
can observe a deviation from the laws of geometrical optics even for the virtual beam. 
Because the incident and transmitted beam are in two different media, the intensity 
is no longer simply proportional to the modulus of the amplitude [28j. To account 
for this distinction we introduce a virtual transmittance V = n\ cost/ cos 9\ : where 
r = sin -1 (sin 0/ri) denotes the angle of transmitted plane waves. As we only consider 
propagating waves in transmission, all the angles and the transmission coefficients are 
real, so that we can ignore the moduli. This also implies that there will be no spatial GH 
shift in this setting. The mean angle for the virtual transmitted beam is thus calculated 
as centroid of the virtual transmitted spectrum: 



JdB 


a(0)\*V(O)r(e) 


Jde 





TV = <r(0)) v = J r ' Ln.L ■ (7) 



Assuming that the spectrum a of the incident beam is symmetric and narrowly 
concentrated around the central plane wave with incidence 9 we expand all 9 dependent 
quantities by setting 9 = 9 + 5 with 5 being small. Identifying t{9q) = r we write 
explicitly 

V{9 + S)^ [1 + S (tan# - r f \ 9o tanr ) + ...], (8) 

cos c/q 

S 2 

r(9 + 5) w r + 5 t'\ 0q + yT"|* + . . . , (9) 

where the prime denotes derivatives with respect to 9. On substitution of these 
expansions into ([7| we can identify the zeroth order in S with Snell's law for the central 
plane wave, while the first correction term is of order 5 2 as the first order terms vanish 
on integration due to symmetry. Collecting all terms of order S 2 yields an approximate 
expression for the mean angle of the virtual beam 



T v ~ T + 



r'|0 o (tan0 o -r\ 0o tanr ) + ^r' 



(5 2 ) = r - V|, o (5 2 ),(10) 
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Figure 3. (Colour online) Construction of the virtual transmitted beam. Every 
incident plane wave is refracted according to Snell's law, which introduces an 
asymmetry in the transmitted spectrum. Averaged over the whole spectrum this 
asymmetry results in an angular shift even for the virtual transmitted beam. The 
figure also shows the coordinate system and angles used in section [3] 



where the equality follows from the properties of r(9) and (S 2 ) is the second spectral 
moment or angular variance of the incident beam calculated as 

2 fd5\a(9 + 5)\W 
{ ° } f d5 \<r(6 + 6)\* ' 1 J 



The second order derivative in ([10]) is not present in other derivations of the angular 
Goos-Hanchen shift on transmission (and it does not occur for reflection [25], [22]). Our 
approach thus differs from known formulas [U [20]. We support our results by the 
excellent agreement with a direct numerical calculation of the virtual shift as defined in 
(see Fig. [i]). The virtual angular GH shift is thus given by 

A« H = TV - T . (12) 

This deviation from Snell's law stems from the asymmetric refraction of the plane waves: 
whereas the incident spectrum is symmetrically distributed around the central plane 
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wave #o, refraction introduces an asymmetry which leads to a shift or deflection of the 
mean angle even for the virtual transmitted beam. 

For the real transmitted beam we have to exchange the virtual transmittance V 
with the real transmittance 7^, a = s(TM), p(TE), which is defined as [27] : 

cosr 



T a (0) = n 



cos 9 



Km 



(13) 



and includes the transmission coefficients. The mean angle of the real reflected beam 
for either s or p polarisation is the centroid of the real transmitted spectrum: 



Ta = (T(0))a = 





a(e)\*T a (e)T(e) 


fde 


°{0)\ 2 T a {0) 



a = s(TM),p(TE). 



In this case we expand the transmittance T a , a = s, p around On 

t' 



T a (0 O + 5) 



COS Th 9 



COS On 



1 + 5 2 



tn 



tan On — 5 t'\$ tan r I + 



(14) 



(15) 



where the prime denotes derivatives with respect to 0. The mean angle of the 
transmitted beam is therefore approximately 



T + 



2/1 



-V 



(s 2 



s(TM),p(TE). 



(16) 



The angular GH shift is the difference between the real and virtual transmitted angle 
A GH = r a — r v for a = s, p, which yields 



AfU) 



f 

AT \n n — 



(5 2 



2cos6» t' 



sin 



2 0ta 



(5 2 



s,p. 



(17) 



This formula is the main result of this section. In this form the angular GH on 
transmission is analogous to its counterpart on reflection ([!]) apart from the geometrical 
factor cos#o/cost . In particular it depends linearly on the second moment of the 
angular spectrum. We also define the total angular GH shift as the sum of (12) and 

(TTrj) as 



A GH + A 



GH 





t'a 

















(5 2 ) 



s,p, 



(18) 



which will be compared to the FF in section [3j 

Of course, the distinction between the angular deflection for the virtual beam and 
the real transmitted beam is to some degree artificial and a measurement of the mean 
angle would give a result corresponding to (18). However, actual measurements are 
often differential, say for different polarisations [8j. In this case the virtual shift cancels 
and the measured quantity corresponds to A GH — A GH . For the purpose of this paper 
it is useful to identify the constituent parts of the angular GH shift. To this end we 
compare in Fig. [i] the total angular GH shift A GH + A GH for both polarisations a = s, p 
with numerical calculations based on (14) for a generic angle of incidence (#o = 0.4rad 
or ^ 23°) of a Gaussian beam with a width w at a glass/air interface (n = 2/3). We 
also show the angular GH shift as defined in (17) and the virtual shift alone. From this 
plot it is obvious that the virtual and real angular GH shift are competing effects, which 
leads to cancellations for the total shift. 
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Figure 4. (Colour online) Angular GH shift at a generic angle of incidence. The plot 
compares the angular GH shift with FF for both s (red) and p (blue) polarisation at 
an incidence angle Oo = 0.4rad (« 23°) and with n = 2/3. The GH shifts are broken 
down into the total shift A™ + A™ (solid), the virtual shift A™ (dashed-dotted) 
and the angular GH shift A™ (dotted). The formulas for the FF shifts A^ F (dashed) 
are derived in section 03 



In [3] the magnitude of the FF at critical incidence is given explicitly in terms of 
the beam waist. In contrast to generic incidence, for the critical angle the effect is not 
inversely proportional to the square of the beam waist, but only to the square root of it. 
It is interesting to see whether the angular GH shift reproduces this result. However, 
because right at the critical angle both the transmission coefficients Q [5]) and the angle 
of refraction ^ change from real to complex functions, the expansions g and (JT5J) 



change in nature. On setting 9 = 9 C — 5 and restricting ourselves to positive angles 
below the critical angle, we can expand the refracted angle in fractional powers of 5 



^> = f-fe^ +0[( ^ )31 ' (19) 

The transmittance can also be expanded in fractional powers of 5 at the critical angle, 



though because T a is present in both the numerator and denominator of (16), and 
because we consider only the first correction term, the coefficients cancel. 



Within the integral (16) the zeroth order of (19) gives rise to Snell's law at critical 
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incidence, namely an angle of refraction of tt/2. The first correction term is essentially 
given by the term of order in the expansion (19), but picks up another from the 
expansion of the transmittance. The mean angle of the transmitted beam is thus given 
by: 

7T r^T r d6 \a(6 c - 5)\ 2 5 , s , s 

^rfa ^iUp «^(TM),p(TE), (20) 

which holds for both s (TM) and p (TE) polarisation. For a general spectrum a 
meaningful interpretation in terms of moments of the spectrum is no longer possible in 
this form, but for a Gaussian beam with spectrum a{8) = exp[—(kwo) 2 (9 c — 5) 2 /4], the 
integrals can be calculated. The angular GH shift for transmission at critical incidence 
is given by 



< A ™ + A ^-SiWi)7t (21) 

where T is the Gamma function [26j . 

Comparing the result with Tureci and Stone [3j (see section [3]) we find the same 
dependence on the square root of the width of the beam, but the proportionality factor is 
different. However, the difference is fairly small, as the factor r ^/4) ~ an< ^ therefore 
close to unity. In their original paper Tureci and Stone state a FF effect of around 30° 
at the critical incidence for a refractive index of n = 1/1.56 and a scaled beam width of 
kw = 8.82. Evaluating (21) for these parameter values gives an angular deflection of 
about 29° which is in line with the different factors. 

This observation is interesting as the difference in the definition of FF and angular 
GH as peak and mean angle is of particular importance at critical incidence as only 
half of the incident spectrum is transmitted. The transmitted spectrum is therefore no 
longer Gaussian, which would enhance the difference between peak and mean angle. A 
fuller account of the shifts at the critical angle which could explain this behaviour is in 
preparation. 

In Fig. [5| we show a comparison between our formula (21) and a numerical 



calculation of (|20j) for two different refractive indices. The agreement is slightly worse 
than for general incidence, but nevertheless excellent. The fact that the shift is negative 
and does not depend on the incident polarisation suggest that the angular GH shift at 
critical incidence is predominantly given by the virtual shift. 



3. Fresnel filtering 

In this section, we study the beam shift defined for the field and not for the intensity as 
in the previous section. The deviation from Snell's law for the field amplitude has been 
studied by Tureci and Stone, who named this effect 'Fresnel filtering' (FF) [3]. To make 
this paper self-contained, we briefly review their theory. We consider a situation where 
a Gaussian beam with central incident angle 9 is propagating from a thicker medium 
with refractive index ri\ to a medium with refractive index 1, such that n = 1/ni, 
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Figure 5. (Colour online) Angular GH shift at a critical incidence. The plot compares 
the angular GH shift with FF at critical incidence C = arcsin(n) for two different 
refractive interfaces: glass - air with n = 2/3 (purple) and GaAs - air with n = 1/3.3 
(orange). The shift is identical for both s and p polarisation. The plot compares the 
total angular GH shift A™ + A™ (solid) with the shift due to FF A£ F (dashed). At 
critical incidence the shift is 3 orders of magnitude larger compared with the shifts in 
Fig.g 



and scattered by an infinite planar interface dividing the two media. We adopt the 
coordinate systems (x^Zi) and {x t) z t ) attached to the incident and refracted beams, as 
shown in Fig.[3| Within the coordinate system Z[) : the incident Gaussian beam with 
a beam waist w is given by 

2 

+ ikz[ 



^i(xi,Zi) = -^rr^exp 



w(zi) 



X 



w(z{) 



(22) 



where w 2 (z\) 



i^-. We assume that the position of the beam waist is located at 



the interface. 

By virtue of the angular spectrum representation, this beam can be expressed as a 
superposition of plane waves 

kw E 



^i{X U Zi) 



2a/^ 



ds a(s) exp [ik {x\ sin 5 + z\ cos S)\ , 



(23) 
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where 5 depends on s via sin 5 = s and cos 5 = y/l — s 2 . The angular spectrum a is 
given by a(s) = exp[— (kw s/2) 2 ]. The transmitted beam is then given by 

^a(x t ,Zt)= 2 ^-° J ds t a (s) a (s)exp [ink (x t sin e + z t cose)] , (24) 

for a = s(TM),p(TE) polarisation. Here, 5 and e are related through sin(# + S) = 
nsin(r + e), which also defines e as a function of s, and i a (s) is the appropriate Fresnel 
transmission coefficient (pi [5]) with 0(s) = 9 + 5(s). On using the polar coodinates 



x t = psin(r — r ) and z t = pcos(r — r ) the refracted beam (24) can be rewritten as 

kiv E f 

ip a (p,r) = 2 ^-° y dst a (s)a(s)exp[ifcpcos(r-r -e))], (25) 

where e depends on s via Snell's law. In the far field, that is in the limit p —> oo, 
the integrand is highly oscillating and can be evaluated using the saddle-point method. 
As pointed out by Tureci and Stone the relevant saddle point is where the cosine in 
the integrand takes its maximum. This relates the free variable r to tq + e and via 
Snell's law also to 9 + 5. For a given r the saddle point s is thus determined by 
solving sin[# + S(so)] = ns'mr for sq- Upon expansion of the exponent and subsequent 
integration the transmitted field is given by 

Up, t) k t a (s )a(s ) J(r; So )e^ (26) 

in the far field. The term J(r; sq) is defined as 



J(r; 5 ) = n 2 . 2 V 1 " 5 o = n ^ V 1 " (27) 

V 1 - n 2 sm^r v cose/ v 

Eq. (26) describes the transmitted field as a function of r and for a saddle point s 
which, for a given r, is fully determined by the incident angle 9q and the refractive 
index n. 



From (26), one can estimate the peak position of the refracted beam for certain 
cases. First, we consider the case for the critical incidence (i.e. 9q = C ), which has 
already been treated by Tureci and Stone [3]. We put r = tt/2 — e with e being a small 
parameter. From sin[# + S(so)] = ncose, we have 

ne 2 

s ~- ^ 7= =■ (28) 



2Vl 



n 



Inserting this into t a (so), &(so), and J(r;so), we get for both p (TE) and s (TM) 
polarisation the same proportionality in the leading order in e 

n 2 ( kwo x 2 



t a (s )cr(so)J{r]So) oc eexp 



e 4 



(29) 



4(1 -n 2 \ 

The different proportionality factor for p and s polarisation does not affect the maximum 
of this expression. Finding the maximum determines the peak position of the electric 
field as 7r/2 — e; the magnitude of the FF is then given by the difference between this 
peak position and the angle expected from Snell's law: 



n 2 J V kw V tan 9 C ^kw 
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This is Tureci and Stone's original result and shows that the angular shift scales as 
(kw^ 1 / 2 at critical incidence. A comparison with the corresponding expression for the 
angular GH shift at critical incidence (21) is shown in Fig. [5| 

For generic incidence Tureci and Stone did not give an explicit analytic result which 
could be compared to the sum A^ H + A^ H [see (12) and (17)]. In the remainder of this 
section we therefore derive an expression for the magnitude of the Fresnel filtering for 
generic incidence. Let us assume that the incident beam is narrow and the incident 
angle is far from the critical angle (i.e. S <C C — #o), so that the profile of the refracted 
beam can be expected to be well approximated by a Gaussian distribution. Assuming 
that the refracted beam is also narrow, we put r = r + e with e being a small parameter. 
Inserting this into sin[# + S(sq)] = nsinr, we have s w (cost / cos9o)ne. Thus, the 
Gaussian distribution a(s ) is written as 

2 



a(so) & exp 



COS T kWQ 

cos 9n 2 



2 2 

n e 



(31) 



Now, let us consider how the functions J(r; so) an d t a (so) shift the above Gaussian 
distribution. To do this, we expand In J(r; s ) an d ln£ a («s ) in terms of e as follows: 



j(0) +j (l) c . 



In J 



Then, we have 



lnt a (s )a(s )J(T;s ) 



(2) 



cos To nkwo 
cos # 2 



(tf>+j (1) )* + 



(32) 
(33) 



(34) 



Ignoring the terms higher than e 2 and completing the square on the right hand side, we 
can rewrite the above equation as 



i » / \ / \ t / \ n f cos t nkw \ ( FF ^2 
\nt a (s )a(s )J(r;s ) « C- — (e - A**) 

\ COS c/q Z J 

where C is a constant and for kw 3> 1. Here, A FF is given by 



A FF = 2 



cos 6 

COS T 



1 



nkwn 



(35) 



(36) 



Eq. (35) tells us that the peak position of the refracted beam is shifted by A FF from 
the transmission angle tq. The coefficients t^ and can be written as 



tn 



n 



s(e=0) 



COST t' a 

COS 9() t n 



5=0 



and 



J 



(i) = i_ 



= -(l-n 2 ) 



S (e=0) 



tan To 
cos 2 6*o 



(37) 



(38) 



where we have changed to the incident angle variable S for later comparison. It is 
interesting to note that both expressions for the angular deflection are a sum of two 
competing effects: tft + in Eq. (36) for FF and A™ + A™ from Eqs. (12) and 
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(17). In fact, on considering all prefactors the term containing is common to both 
formulas. Nevertheless the results can be quite different as shown in Fig. [4j where the 
formulas for FF and the angular GH shift predict deflections in opposite directions. 

4. Discussion 

In this paper we have highlighted the similarities and differences between Fresnel filtering 
(FF) and the angular Goos-Hanchen shift (GH) on transmission: two related effects 
which give rise to a deviation from Snell's law if the incident field is not a plane wave, 
but a spatially localised beam of light. Both effects are based on a modification or 
'filtering' of the angular spectrum, and are therefore of common in nature. However, 
different assumptions in the definition of the two effects gives rise to a difference in 
the magnitude of the deflections. To examine these differences we have extended the 
existing theories for FF and the angular GH shift. For the former we have derived an 
analytical expression for the deviation at generic incidence, that is not for the special 
case of critical incidence, and for the latter we have found an approximate analytic result 
for the special case of critical incidence. 

Before we discuss the results in detail we clearly state the differences in definition. 
FF has been defined as the difference between the peak transmission angle in far field 
and the angle expected from Snell's law, whereas the angular GH shift on transmission 
compares Snell's law with the mean angle. As a consequence FF uses the amplitude 
transmission coefficient, whereas the angular GH shift depends on the transmittance, 
that is transmission coefficient for the energy or intensity. 

In light of these differences it may seem surprising to find any similarities, but 
of course amplitude and intensity, near field and far field are related and we do find 
terms common to both formulas. Most prominently perhaps the logarithmic derivative 



of the transmission coefficient in Eqs. (17) and (36), which is a familiar term from 



reflection. Interestingly, it is not this shared term which explain the very similar result 
for critical incidence. At this special angle the deviation from Snell's law calculated in 
both theories is almost identical, and the remaining difference can be explained as the 



difference between mean and peak angle. In terms of the competing terms in (18) and 



(36), however, it is not the common term which is dominant at the critical angle, but, 
somewhat surprisingly, the 'virtual' shift A^ H or the term. 

For generic incidence the magnitudes of thrangular GH shift and FF can be quite 
different, though of course the precise difference depends on the angle of incidence. Our 
choice of 9 showcases the differences. Angles of incidence closer to critical incidence 
lead to smaller differences. As the logarithmic derivative of the transmission coefficient 
is common to both, it is again the 'virtual' term which is responsible for the difference. 
In case of the FF this term is always negative (for n < 1) and larger in magnitude which 
explains the overall negative shift. For the angular GH shift the virtual term is also 
always negative if n < 1, but smaller than the derivative term, which is why the shift 
turns negative for larger angles of incidence. The fundamental difference may be founded 
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in the distinction between amplitude and field, which raises the question whether the 
difference between angular GH shift and FF has it roots in the difference between the 
electric field and the displacement field. An answer to this question, however, is outside 
the scope of this paper and deserves its own exposition. 
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